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Abstract 

We study the short-range nucleon-nucleon interaction in a chiral constituent 
quark model by diagonalizing a Hamiltonian comprising a linear confinement 
and a Goldstone boson exchange interaction between quarks. The six-quark 
harmonic oscillator basis contains up to two excitation quanta. We show that 
the highly dominant configuration is | s 4 p 2 [42]o[51]fs > due to its specific 
flavour-spin symmetry. Using the Born-Oppenheimer approximation we find a 
strong effective repulsion at zero separation between nucleons in both 3 Si and 
1 So channels. The symmetry structure of the highly dominant configuration 
implies the existence of a node in the S-wave relative motion wave function at 
short distances. The amplitude of the oscillation of the wave function at short 
range will be however strongly suppressed. We discuss the mechanism leading 
to the effective short-range repulsion within the chiral constituent quark model 
as compared to that related with the one-gluon exchange interaction. 



I. INTRODUCTION 



An interest in the constituent quark model (CQM) has recently been revitalized Q after 
recognizing the fact that the constituent (dynamical) mass of the light quarks appears as a 
direct consequence of the spontaneous chiral symmetry breaking (SCSB) and is related 
with the light quark condensates < qq > of the QCD vacuum. This feature becomes explicit 
in any microscopical approach to SCSB in QCD, e.g. in the instanton gas (liquid) model [|J. 
The mechanism of the dynamical mass generation in the Nambu-Goldstone mode of chiral 
symmetry is very transparent within the a-model |J or Nambu and Jona-Lasinio model 
||. Another consequence of the chiral symmetry in the Nambu-Goldstone mode is the 
appearance of an octet of Goldstone bosons (n,K,r] mesons). It was suggested in |jj that 
beyond the scale of SCSB, nonstrange and strange baryons should be viewed as systems 
of three constituent quarks which interact via the exchange of Goldstone bosons and are 
subject to confinement. This type of interaction between the constituent quarks provides a 
very satisfactory description of the low-lying nonstrange and strange baryon spectra |l],|7||i 
including the correct ordering of the levels with positive and negative parity in all parts of 
the considered spectrum. 

So far, all studies of the short-range NN interaction within the constituent quark model 
were based on the one-gluon exchange interaction (OGE) between quarks. They explained 
the short-range repulsion in the NN system as due to the colour-magnetic part of OGE 
combined with quark interchanges between 3q clusters. (For reviews and earlier references 
see [PHTT||). There are also models which attribute the short-range repulsion in the NN 



system to the colour-electric part of OGE [12 



In order to provide the necessary long- and intermediate-range attraction in the baryon- 
baryon system, hybrid models were suggested fl3|-|T5f , where in addition to OGE, the quarks 
belonging to different 3q clusters interact via pseudoscalar and scalar meson exchange. In 
these hybrid models the short-range repulsion in the NN system is still attributed to OGE 
between the constituent quarks. 
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It has been shown, however |p],^6|, that the hyperfine splittings as well as the correct 
ordering of positive and negative parity states in spectra of baryons with u,d,s quarks are 
produced in fact not by the colour-magnetic part of OGE, but by the short-range part of 
the Goldstone boson exchange (GBE) interaction. This short-range part of GBE has just 
opposite sign as compared to the Yukawa tail of the GBE interaction and is much stronger 
at short interquark separations. There is practically no room for OGE in light baryon 
spectroscopy and any appreciable amount of colour-magnetic interaction, in addition to 
GBE, destroys the spectrum |§. If so, the question arises which interquark interaction is 
responsible for the short-range NN repulsion. The goal of this paper is to show that the same 
short-range part of GBE, which produces good baryon spectra, also induces a short-range 
repulsion in the NN system. 

The present study is rather exploratory. We calculate an effective NN interaction at 
zero separation distance only. We also want to stress that all main ingredients of the NN 
interaction, such as the long- and middle-range attraction and the short-range repulsion 
are implied by the chiral constituent quark model. Indeed, the long- and middle-range 
attraction automatically appear in the present framework due to the long-range Yukawa tail 
of the pion-exchange interaction between quarks belonging to different nucleons and due to 
2tt (or sigma) exchanges. Thus, the only important open question is whether or not the 
chiral constituent quark model is able to produce a short-range repulsion in the NN system. 

For this purpose, we diagonalize the Hamiltonian of Ref. in a six-quark harmonic 
oscillator basis up to two excitations quanta. Using the Born-Oppenheimer (adiabatic) 
approximation, we obtain an effective internucleon potential at zero separation between 
nucleons from the difference between the lowest eigenvalue and two times the nucleon mass 
calculated in the same model. We find a strong effective repulsion between nucleons in 
both 3 Si and 1 Sq channels of a height of 800-1300 MeV. This repulsion implies a strong 
suppression of the NN wave function in the nucleon overlap region as compared to the wave 
function of the well separated nucleons. 

Due to the specific flavour-spin symmetry of the GBE interaction, we also find 
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that the highly dominant 6q configuration at zero-separation between nucleons is 
|s 4 p 2 [42] [51]i?5 >. As a consequence the 6q region (i.e. the nucleon overlap region) can- 
not be adequately represented by the one-channel resonating group method (RGM) ansatz 
A{N(1, 2, 3)iV(4, 5, 6)%(r)} which is commonly used at present for the short-range NN in- 
teraction with the OGE interaction. 

The symmetry structure [42]o[51]_fs of the lowest configuration will induce an additional 
effective repulsion at short range related to the "Pauli forbidden state" in this case. This 
latter effective repulsion is not related to the energy of the lowest configuration as com- 
pared to two-nucleon threshold and thus cannot be obtained within the Born-Oppenheimer 
approximation procedure. We notice, however, that the structure of the six-quark wave 
function in the nucleon overlap region is very different from the one associated with the soft 
or hard core NN potentials. 

This paper is organized as follows. In section 2, in a qualitative analysis at the Casimir 
operator level, we show that the short-range GBE interaction generates a repulsion between 
nucleons in both 3 Si and l So channels. We also suggest there that the configuration with 
the [51]^ flavour-spin symmetry should be the dominant one. Section 3 describes the 
Hamiltonian. Section 4 contains results of the diagonalization of the 6q Hamiltonian and 
of the NN effective interaction at zero separation between nucleons. The structure of the 
short-range wave function is also discussed in this section. In section 5, we show why the 
single-channel RGM ansatz is not adequate in the present case. In section 6, we present a 
summary of our study. 

II. A QUALITATIVE ANALYSIS AT THE CASIMIR OPERATOR LEVEL 

In order to have a preliminary qualitative insight it is convenient first to consider a 
schematic model which neglects the radial dependence of the GBE interaction. In this model 
the short-range part of the GBE interaction between the constituent quarks is approximated 
by i 
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V x = -C x J2\f.Xfa i .a j , (1) 

i<j 

where X F with an implied summation over F (F=l,2,...,8) and a are the quark flavour Gell- 
Mann and spin matrices respectively. The minus sign of the interaction ([!]) is related to the 
sign of the short-range part of the pseudoscalar meson-exchange interaction (which is oppo- 
site to that of the Yukawa tail), crucial for the hyperfine splittings in baryon spectroscopy. 
The constant C x can be determined from the A — N splitting. For that purpose one only 
needs the spin (S), flavour (F) and flavour-spin (FS) symmetries of the N and A states, 
identified by the corresponding partitions [f] associated with the groups SU(2) S , SU(3)p 
and SU(6)fs '■ 

\N> = \s 3 [3} FS [21] F [21} s >, (2) 
|A > = |s 3 [3] FS [3] F [3] 5 > . (3) 

Then the matrix elements of the interaction (H) are Jj]] : 

< N\V X \N > = -UC X , (4) 

< A\V X \A > = -4C X . (5) 

Hence — En = 10C X , which gives C x = 29.3MeV, if one uses the experimental value of 
293 MeV for the A — iV splitting. 

To see the effect of the interaction (JI|) in the six-quark system, the most convenient is to 
use the coupling scheme called FS, where the spatial [f]o and colour [f] c parts are coupled 
together to [/]oc> an d then to the SU(6)fs flavour-spin part of the wave function in order to 
provide a totally antisymmetric wave function in the OCFS space |i~8 |. The antisymmetry 
condition requires [/]fs = [/]oc> where [/] is the conjugate of [/]. 

The colour-singlet 6q state is [222] c- Assuming that N has a [3]o spatial symmetry, 
there are two possible states [6]q and [42] q compatible with the S-wave relative motion in 
the NN system [[L7|]. The flavour and spin symmetries are [42] ^ and [33] 5 for 1 S and [33] .p 
and [42] 5 for 3 Si channels. Applying the inner product rules of the symmetric group for 



both the [f]o x [f]c and [f] F x [/]§ products one arrives at the following 6q antisymmetric 
states associated with the 3 Si and 1 5'o channels []T8| , [T9[1 : |[6]o[33]^,g >, | [42] [33]i7s >, 
|[42] [51] FS >, |[42]o[411] FS >, |[42] [321] FS >, |[42] [2211] FS >. 

Then the expectation values of the GBE interaction ([I]) for these states can be easily 
calculated in terms of the Casimir operators eigenvalues for the groups SU (6) ps, SU(3) F and 
SU(2)s using the formula given in Appendix A. The corresponding matrix elements are given 
in Table 1, from where one can see that, energetically, the most favourable configuration is 
[51]fs- This is a direct consequence of the general rule that at short range and with fixed spin 
and flavour, the more "symmetric" a given FS Young diagram is, the more negative is the 
expectation value of ([!]) |IJ. The difference in the potential energy between the configuration 
[51]irg and [33] .ps or [411]p 5 is of the order : 

< [33] FS \V X \[33] FS > - < [51] FS \V X \[51] FS > = 

(6) 

< [411] FS |Vy[411] FS > - < [5l} FS \V x \[5l} FS > = 2AC X 
and using C x given above one obtains approximately 703 MeV for both the SI = 10 and 01 
sectors. 

In a harmonic oscillator basis containing up to 2hu excitation quanta, there are two 
different 6q states corresponding to the [6]o spatial symmetry with removed center of mass 
motion. One of them, |s 6 [6]o >, belongs to the N = shell, where N is the number of 
excitation quanta in the system, and the other, ^\s 5 2s[6]o > — y|l s4 P 2 [6]o >, belongs 
to the N = 2 shell. There is only one state with [42] o symmetry, the |s 4 p 2 [42]o > state 
belonging to the N=2 shell. While here and below we use notations of the shell model it is 
always assumed that the center of mass motion is removed. 

The kinetic energy KE for the |s 4 p 2 [42]o > state is larger than the one for the |s 6 [6]o > 
state by KEn=2 — KEn=o = hu. Taking Tiuj ~ 250 MeV and denoting the kinetic energy 
operator by H Q , we obtain : 

< s 6 [33} FS \H + V x \s 6 [33] FS > - < s 4 p 2 [51] FS \H + V x \sY[51] FS >~ 453MeV (7) 

which shows that [51]^ is far below the other states of Table 1. For simplicity, here we 
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have neglected a small difference in the confinement potential energy between the above 
configurations. 

This qualitative analysis suggests that in a more quantitative study, where the radial 
dependence of the GBE interaction is taken into account, the state |s 4 p 2 [42] [51]fs > will 
be highly dominant and, due to the important lowering of this state by the GBE interaction 
with respect to the other states, the mixing angles with these states will be small. That this 
is indeed the case, it will be proved in the section 4 below. 

Table 1 and the discussion above indicate that the following configurations should be 
taken into account for the diagonalization of the realistic Hamiltonian in section 4: 

|1 > = \s 6 [6]o[S3] FS > 

|2 > = |sV [42] G [33] FS > 

|3 > = |sV[42]o[51]fs > 

( 8 ) 

|4> = |sV[42]o[411]fs > 

|5> = |(yf S 5 2 S -yi S V)[6]o[33] FS > 

A strong dominance of the configuration |3 > also implies that the one-channel approxi- 
mation A{NNx(r)} is highly inadequate for the short-range NN system. This problem will 
be discussed in section 5. 

Now we want to give a rough estimate of the interaction potential of the NN system 
at zero separation distance between nucleons. We calculate this potential in the Born- 
Oppenheimer (or adiabatic) approximation defined as : 

V NN (R) =<H> R -<H >OQ (9) 

where R is a collective coordinate which is the separation distance between the two s 3 
nucleons, < H >r is the lowest expectation value of the Hamiltonian describing the 6q 
system at fixed R and < H >oo= 2m N for the NN problem, i.e. the energy of two well 
separated nucleons. As above we ignore the small difference between the confinement energy 
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of < H >r=o and < H > OQ . That this difference is small follows from the A?.A^ structure of 
the confining interaction and from the identity : 

< [222] c |£ A l c .A||[222] c >= 2 < [lll] c | £ Af-A^lll], > . (10) 

If the space parts [6]o and [3]o contain the same single particle state, for example an s-state, 
then the difference is identically zero. 

It has been shown by Harvey |TB[ that when the separation R between two s 3 nucleons 
approaches 0, then only two types of 6q configurations survive: |s 6 [6]o > and |s 4 p 2 [42]o >• 
Thus in order to extract an effective NN potential at zero separation between nucleons 
in adiabatic Born-Oppenheimer approximation one has to diagonalize the Hamiltonian in 
the basis |1 > — 14 >. In actual calculations in section 4 we extend the basis adding 
the configuration |5 >, which practically does not change much the result. For the rough 
estimate below we take only the lowest configuration |3 >. One then obtains 

< s A p 2 [42]o[51]fs\H + V x \s 4 P 2 [42]o[51]fs > -2 < N\H + V X \N >= 

-100/3 + 28)C X + 7/Ahcu = 280 MeV, if SI = 10 
(-32 + 28)C X + 7/Ahcu = 320 MeV, if SI = 01 

The rough estimate (O) suggests that there is an effective repulsion of approximately equal 
magnitude in the NN system in the nucleon overlap region in both 3 Si and 1 S'o channels. In 
a more quantitative calculation in Section 4 we find that the height of the effective core is 
much larger, in particular 830 MeV for 3 Si and about 1.3 GeV for 1 S'o. 

At this stage it is useful to compare the nature of the short-range repulsion generated 
by the GBE interaction to that produced by the OGE interaction. 

In the constituent quark models based on OGE the situation is more complex. Table 1 
helps in summarizing the situation there. In this table we also give the expectation value of 
the simplified chromo-magnetic interaction 

V cm = -C cm ^2 (12) 



in units of the constant C cm (the constant C cm can also be determined from the A — N 
splitting to be C cm ~ 293/16 MeV). 

The expectation values of ( |i~2"|) can be easily obtained in the CS scheme with the help 
of Casimir operator formula in Appendix A and can be transformed to FS scheme by using 
the unitary transformations from the CS scheme to the FS scheme given in Appendix B. 

The colour-magnetic interaction pulls the configuration |s 4 p 2 [42]o[42]cs > down to be- 
come approximately degenerate with |s 6 [6]o[222]cs > which is pulled up. In a more detailed 
calculation with explicit radial dependence of the colour-magnetic interaction as well as with 
a Coulomb term the configuration |s 6 [6]o > is still the lowest one |],[2(J. (With the model 
(12]) the hu should be about 500 MeV). Thus in the Born-Oppenheimer approximation we 
can roughly estimate an effective interaction with OGE model through the difference 

< s 6 [6} o [222} cs \H + V cm \s 6 [6} [222]cs > -2 < N\H + V cm \N > 



C cm + 3/4hu = IVlMeV if SI = 10 



3 (13) 
24C cm + 3/4hu = 815MeV if SI = 01 

We conclude that both the GBE and OGE models imply effective repulsion at short 
range of approximately same magnitude. 

III. THE HAMILTON! AN 



In this section we present the GBE model [[y,[7| used in the diagonalization of six-quark 
Hamiltonian in the basis (|8|). The Hamiltonian reads 



H = 6m + Y / ^-- + E W'y) + E (14) 

i i<j Kj 

where m is the constituent quark mass; = |f, — f}| is the interquark distance. 
The confining interaction is 

V conf (r tj ) = -l^-X c 3 Cr tJ (15) 



where A? are the SU(3)-colour matrices and C is a parameter given below. 

The spin-spin component of the GBE interaction between the constituent quarks i and 
j reads: 



mo 



(iX(^)AfAj 
[f=i 



+ E Vk(fy)Af Af + V^(^-)AfA| + K/(^)A°A° * ■ (16) 

where A F , F = 1, 8 are flavour Gell-Mann matrices and A = ^2/3 1, where 1 is the 3x3 
unit matrix. Thus the interaction ( |T6"D includes 7T, K, rj and r/ exchanges. While the it, K, 
r) mesons are (pseudo)Goldstone bosons of the spontaneously broken SU(3)l x SU(3)r — > 
S77(3)y chiral symmetry, the 7/ (flavour singlet) is a priori not a Goldstone boson due to the 
axial U(1)a anomaly. In the large Nq limit the axial anomaly disappears, however, and the 
rj' becomes the ninth Goldstone boson of the spontaneously broken U{3)l x U(3)r — ► t/(3)y 
chiral symmetry |^T|. Thus in the real world with Nq = 3 the rj' should also be taken 
into account, but with parameters essentially different from tt, K, rj exchanges due to 1/Nq 
corrections. For the system of u and d quarks only the .fT-exchange does not contribute. 

In the simplest case, when both the constituent quarks and mesons are point-like particles 
and the boson field satisfies the linear Klein-Gordon equation, one has the following spatial 
dependence for the meson-exchange potentials M : 

a 2 1 1 e -H r ij 

^■) = g3^W— (7 = ( 17 ) 

where /x 7 are the meson masses and g^/4:7i are the quark- meson coupling constants given 
below. 

Eq. (p~7| ) contains both the traditional long-range Yukawa potential as well as a 5-function 
term. It is the latter that is of crucial importance for baryon spectroscopy and short-range 
iVTV" interaction since it has a proper sign to provide the correct hyperfine splittings in 
baryons and is becoming highly dominant at short range. Since one deals with structured 
particles (both the constituent quarks and pseudoscalar mesons) of finite extension, one must 
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smear out the 5-function in flPT|). In Ref. Jj]] a smooth Gaussian term has been employed 
instead of the 5-function 

4 

4m8{fij) => -^a 3 exp(-ct 2 (r - r ) 2 ). (18) 

'7T 



where a and r are adjustable parameters. 

The parameters of the Hamiltonian fll4|) are pj: 



2 2 2 

%i = %i = 0.67; ^ = 1.206 
An Aix Aix 

r = 0.43 fm, a = 2.91 far 1 , C = 0.474 fai" 2 , m = 340 MeV. 

^ = 139 MeV, fx v = 547 MeV, /v = 958 MeV. (19) 

The Hamiltonian ([14]) with the parameters fll9|) provides a very satisfactory description of 
the low- lying A" and A spectra in a fully dynamical nonrelativistic 3-body calculation [^]. 

At present we are limited to use a |s 3 > harmonic oscillator wave function for the nucleon 
in the NN problem. The parametrization (|I~9"D is especially convenient for this purpose since 
it allows to use the |s 3 > as a variational ansatz. Otherwise the structure of should be 
more complicated. Indeed, < A^|if|A^ > takes a minimal value of 969.6 MeV at a harmonic 



oscillator parameter value of (3 = 0.437 fm |[22|| , i.e. only 30 MeV above the actual value 



in the dynamical 3-body calculation. In this way one satisfies one of the most important 
constraint for the microscopical study of the A^V interaction : the nucleon stability condition 
1 

d 

— <N\H\N>=0. (20) 
do 



The other condition, the qualitatively correct A — N splitting, is also satisfied j22|. 

We keep in mind, however, that a nonrelativistic description of baryons cannot be com- 
pletely adequate. Within the semirelativistic description of baryons |J the parameters 
extracted from the fit to baryon masses become considerably different and even the repre- 
sentation of the short-range part of GBE (|i~8|) has a different form. Within a semirelativistic 
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description the simple s 3 wave function for the nucleon is not adequate anymore. All this 
suggests that the description of the nucleon based on the parameters fll^D and an s 3 wave 
function is only effective. Since in this paper we study only qualitative effects, related to 
the spin-flavour structure and sign of the short-range part of GBE interaction, we consider 
the present nonrelativistic parametrization as a reasonable framework. 

We diagonalize the Hamiltonian fll4|) in the basis All the necessary matrix elements 
are calculated with the help of the fractional parentage technique. Some important details 
can be found in Appendices C and D. 

IV. RESULTS AND DISCUSSION 

In Tables II and III we present our results obtained from the diagonalization of the 
Hamiltonian (|I4~D in the basis (|J). According to the definition of the effective potential within 
the Born-Oppenheimer approximation (|9|) at zero separation between nucleons all energies 
presented in the Tables II and III are given relative to two-nucleon threshold, i.e. the quantity 
2 < N\H\N >= 1939 MeV has always been subtracted. In the second column we present 
the diagonal matrix elements for all the states listed in the first column. In the third column 
we present all the eigenvalues obtained from the diagonalization of a 5 x 5 matrix. In the 
fourth column the amplitudes of all components of the ground state are given. In agreement 
with Sec. 2, one can see that the expectation value of the configuration |s 4 p 2 [42]o[51]fs > 
given in column 2 is much lower than all the other ones, and in particular it is about 1.5 GeV 
below the expectation value of the configuration |s 6 [6]o[33]fs >• The substantial lowering of 
the configuration |s 4 p 2 [42]o[51]i?5 > relative to the other ones implies that this configuration 
is by far the most important component in the ground state eigenvector. The last column 
shows that the probability of this configuration is 93% both for SI = 10 and SI = 01. As a 
consequence, the lowest eigenvalue is only about 100 MeV lower than the expectation value 
of the configuration |3 >. 

The main outcome is that Vnn(R = 0) is highly repulsive in both 3 Si and 1 S channels, 
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the height being 0.830 GeV in the former case and 1.356 GeV in the latter one. 

In order to see that it is the GBE interaction which is responsible for the short range re- 
pulsion, it is very instructive to remove V x from the Hamiltonian (|HD , compute the "nucleon 
mass" in this case, which turns out to be itln = 1.633 GeV at the harmonic oscillator pa- 
rameter (3 = 0.917 fm and diagonalize such a Hamiltonian again in the basis (^j). In this case 
the most important configuration is |s 6 [6]o[33]i?s >• Subtracting from the lowest eigenvalue 
the "two-nucleon energy" 2m N = 2 x 1.633 GeV one obtains GBE {R = 0) = -0.197 
GeV. This soft attraction comes from the unphysical colour Van der Waals forces related 
to the pairwise confinement. The Van der Waals forces would not appear if the basis was 
restricted to the |s 6 > state only. If the spatially excited 3q clusters from the s 4 p 2 config- 
urations were removed the Van der Waals forces would disappear and we would arrive at 
Vnn GBE {R = 0) = 0. Thus it is the GBE interaction which brings about 1 GeV repulsion, 
consistent with the previous discussion. 

The effective repulsion obtained above implies a strong suppression of the L = relative 
motion wave function in the nucleon overlap region, as compared to the wave function of 
two well separated nucleons. 

There is another important mechanism producing additional effective repulsion in the 
NN system, which is related to the symmetry structure of the lowest configuration but 
not related to its energy relative to the NN threshold. This "extra" repulsion, related 
to the "Pauli forbidden state" 23]] , persists if any of the configurations from the \s A p 2 > 



shell becomes highly dominant JT7|]. Indeed, the NN phase shift calculated with a pure 
[51]^^ state, which is projected "by hands" (not dynamically) from the full NN state in 
a toy model ||, shows a behaviour typical for repulsive potentials. As a result the S-wave 



NN relative motion wave function has an almost energy independent node |24]]. A similar 
situation occurs in 4 He — 4 He scattering ||25|| . The only difference between this nuclear case 



and the NN system is that while in the former case a configuration s 8 is indeed forbidden 
by the Pauli principle in eight-nucleon system, the configuration s 6 is allowed in a six-quark 
system, but is highly suppressed by dynamics, as it was discussed above. In the OGE model 
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this effect is absent because none of the [42] o states is dominant P, 20 ,^7 1 . The existence of 



a strong effective repulsion, related to the energy balance in the adiabatic approximation, 
as in our case, suggests, however, that the amplitude of the oscillating NN wave function at 
short distance will be strongly suppressed. 

To illustrate this discussion we project the lowest eigenvector in Table II onto the NN 
and AA channels. The projection onto any baryon-baryon channel B1B2 is defined as follows 



= < ^(1,2,3)^(4, 5, 6)|¥(1, 2, ..,6) >, (21) 

where 2, ...,6) is a fully antisymmetric 6q wave function, which in the present case is 
represented by the eigenvector in Table II, and Bi(l, 2, 3) and -82(4, 5, 6) are intrinsic baryon 
wave functions. 

In order to calculate ([H]) we need a "3 + 3" expansion of each state in the basis (^). The 



corresponding "3 + 3" decomposition of each state can be found in [2^] in the CS coupling 
scheme. To use it here one needs the unitary matrix from the CS basis to the FS one. This 
matrix can be found in Appendix B. 



In Fig. 1 we show the projections fl2"T|) onto the NN and AA channels in the 3 Si partial 



NN wave at short range. In fact, such projections can be shown for other channels too as 
e.g. NN*, N*N*,... some of them being not small. Note that our six-quark wave function, 
calculated at short range only, was normalized to 1. Hence, we cannot show the suppression 
of the NN projection in the nucleon overlap region as compared to the wave function of 
the well separated nucleons, discussed above. This can only be seen from ^ NN (r) obtained 
in dynamical calculations including not only the short-range 6q configurations, like in the 
present paper, but also the basis states representing the middle- and long distances in the 
NN system. 

In Fig. 1 one observes a nodal behaviour of both ^NN(r) and \1/aa(^) at short range. 
Also \I/aa(^*) is essentially larger at short range than ^/ NN (r). At large distances only ^nn^ 
will survive. This nodal behaviour is related to the fact that the configuration |3 > is highly 
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dominant. In the case of any configuration s 4 p 2 or s 5 2s from the N=2 shell, the relative 
motion of two s 3 clusters (e.g. NN and AA) is described by a nodal wave function. 

Now we want to discuss the question which type of NN potential would be equivalent 
to the short-range picture obtained above. If one considers the effect of the short-range 
dynamics on the NN phase shifts, in a limited energy interval the phase shifts in both 3 Si and 
1 Sq partial waves can be simulated by strong repulsive core potentials or by "deep attractive 
potentials with forbidden states" W% - The latter potentials are in fact supersymmetric 
partners of the former ones 



If, on the other hand, one considers the effect of the short-range dynamics on the structure 
of the wave function at short range, it is difficult to construct a potential which would be 
adequate. For example, a repulsive core potential produces a wave function which is indeed 
suppressed at short range, but does not have any nodal structure. If one takes, instead, 
a "deep attractive potential with forbidden state" , one obtains a nodal behaviour, but the 
wave function is not suppressed at short range (i.e. the amplitude left to the node is a 
very large one). As a direct consequence, the latter potential produces a very rich high- 
momentum component, which is in contradiction with the deuteron electromagnetic form 



factors [2^| . The high-momentum component, implied by a "very soft node", like in our 
case, will be much smaller and closer to that one obtained from the potentials with strong 
repulsive core. 

We also see large projections onto other BiB 2 channels (exemplified by the AA channel 
in Fig. 1). These components cannot be taken into account in any simple NN potential, in 
principle. Thus, if we are interested in effects, related to the short-range NN system, there 
is no way, other than to consider the full 6q wave function in this region. 
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V. WHY THE SINGLE CHANNEL RESONATING GROUP METHOD ANSATZ 

IS NOT ADEQUATE ? 



In this section we show that the currently used one-channel resonating group method 
(RGM) ansatz for the two-nucleon wave function is not adequate in a study of the short- 
range NN interaction with the chiral quark model. 

In the one-channel RGM approximation the 6q wave function has the form 

</> = A{N(1, 2, 3)JV(4, 5, 6) X (f)}, (22) 
A=^=(1-9P 36 ), 



T\ + r 2 + r 3 r 4 + r 5 + r 6 



This is reasonable in the case where the short-range quark dynamics is described in terms 



of the OGE interaction. In this case the addition of new channels, orthogonal to (|22|) , does 
not change considerably the full wave function in the nucleon overlap region. This is not the 
case for the chiral constituent quark model, where the short-range quark dynamics is due 
to the GBE interaction. To have a better insight why (|22|) is a poor approximation in the 
present case, we begin with the explanation why (|22"D is reasonable for OGE model |2C . 



To this end, it is very convenient to use the six-quark shell model basis for the NN function 
in the nucleon overlap region []T75] , pU] , pD|] . Such a basis is much more flexible than 



Diagonalizing a Hamiltonian comprising OGE and a confining interaction in the harmonic 
oscillator basis up to two excitation quanta, one can obtain the 6q wave function in the form 
2| 



$ = Co\s 6 > +J2C a \a >, (23) 

a 

where a lists all possible configurations in the N=2 shell : [6]o[222] cs , [42]o[42]cs, 
[42]o[321] os , [42]o[222] cs , [42] [3111]os, [42] [21111] os . With the OGE interaction, the 
CS coupling scheme based on the chain SU(6)cs ^ SU(3)o x SU(2)s is more convenient. 
It has been found that there are a few most important configurations - |s 6 [6]o[222]cs > 
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, |sV[42] [42] cs >, |sV[42]o[321]c 5 >, \(^5/6s 5 2s - <Jl/6s 4 p>)[G]o[222]cs > - with size- 
able amplitudes C a 

Now, let us expand the RGM wave function (^) in the shell model basis. For that 
purpose, the trial function x(r) in (|22l) should be expanded in a harmonic oscillator basis 
too 

XL=o(r) = < Xl=o\<Pns > 4>Ns(r), (24) 

N =0,2,4,... 

where <j)Ns(r) is a harmonic oscillator state with N quanta and L = 0. Thus in the ansatz 



(p2|) the variational coefficients based on the expansion fl24|) are < Xl=o\4>ns >• The last 
step is to use the expressions (21) and (22) of Ref. @ for A{N(1, 2, 3)JV(4, 5, 6)0 o s(r)}, 
A{iV(l, 2, 3)iV(4, 5, 6)025(0}, written in the shell model basis. These are transformations 
from one basis to another and do not depend on the 6q dynamics. If it turns out that for 
a given Hamiltonian the variational coefficients C a in (^) are close to the algebraical ones 
< a\A{N(l, 2, 3)iV(4, 5, 6)(/>2s(r)} >, then one can conclude that P2] ) is a good approxima- 
tion for the variational solution (1231) . If not, the variational ansatz (EI3) is poor and other 



channels, not equivalent to (H), should be added (e.g. A{N*Nx*(r)},A{N*N*x**(r)},...). 



For the OGE model it is found that indeed the variational coefficients C a in (|23D are very 
close to the algebraical ones |2(| (see also HU). 

Let us now turn to the analysis of the results of Sec. 4 based on the GBE interaction. 
Using the unitary transformation from the CS to FS scheme, given in Appendix B, one can 



rewrite Eqs.(21) and (22) of Ref. |26 



as : 



i{iV(l,2,3)iV(4,5,6)0os(r)} 5 / = io = J y |s 6 [6] [33] FS >, (25) 

i{iV(l,2,3)iV(4,5,6)0 2s (f)} s/=1 o = ^\(^s 5 2s - y| S V)[6] [33] FS > 

Ay/2, 4 2Nol rool ^ 4V2, 4 2r 



o ,sV[42] [33] FS > ^-|sV[42]o[51]f5>- (26) 

From the expression (pH|) we see that the relative amplitudes of the states |5 >, |2 > and 
1 3 > are in the ratio 
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|5 >: |2 >: |3 >= 3 : -4 : -4 (27) 

and the amplitude of the state |4 > is zero. The diagonalization of the Hamiltonian made 
in the previous section gives 

|5 >: \2 >: |3 >: |4 >~ 0.06 : 0.08 : -0.96 : 0.20, (28) 

Therefore the ansatz (^) is completely inadequate in the nucleon overlap region and the 
incorporation of additional channels is required in RGM calculations. 

VI. SUMMARY 

In the present paper we have calculated an adiabatic NN potential at zero separation 
between nucleons in the framework of a chiral constituent quark model, where the constituent 
quarks interact via pseudoscalar meson exchange. Diagonalizing a Hamiltonian in a basis 
consisting of the most important 6q configurations in the nucleon overlap region, we have 
found a very strong effective repulsion of the order of 1 GeV in both 3 Si and 1 Sq NN 
partial waves. Due to the specific flavour-spin symmetry of the Goldstone boson exchange 
interaction the configuration |s 4 p 2 [42]o[51].F5 > becomes highly dominant at short range. As 
a consequence, the projection of the full 6q wave function onto the NN channel should have 
a node at short range in both 3 Si and 1 5'o partial waves. The amplitude of the oscillation 
left to the node should be strongly suppressed as compared to the wave function of two well 
separated nucleons. 

We have also found that due to the strong dominance of the configuration 
|s 4 p 2 [42]o[51]i?5 > the commonly used one-channel RGM ansatz is a very poor approxi- 
mation to the 6q wave function in the nucleon overlap region. 

Thus, within the chiral constituent quark model one has all the necessary ingredients to 
understand microscopically the NN interaction. There appears strong effective short-range 
repulsion from the same part of Goldstone boson exchange which also produces hyperfine 
splittings in baryon spectroscopy. The long- and middle-range attraction in the NN system 
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is automatically implied by the Yukawa part of pion exchange and two-pion (or a) exchanges 
between quarks belonging to different nucleons. With this first encouraging result, it might 
be worthwhile to perform a more elaborate calculation of NN system and other baryon- 
baryon systems within the present framework. 

APPENDIX A: 

The expectation value of the operators ([[]) and (|I2]), displayed in Table [I], are calculated 
with the following formulae: 

< £ Xi.Xj&i.&j >= AC S 2 U (6) - 2Cf /(3) - Icf {2) - 8N (Al) 

i<j 3 

where N is the number of particles, here N=6, and C 2 U ^ is the Casimir operator eigenvalues 
of SU(n) which can be derived from the expression : 

C S 2 U(n) = \[f[{f[ + n - 1) + f 2 {f 2 + n - 3) + &(f 3 + n - 5) 

i n— 1 

+Wi + n - 7) + ... + - n + 3)] - — (E (A2) 

where /• = fi — f n , for an irreductible representation given by the partition [/i, / 2) /n]- 

APPENDIX B: 

This appendix reproduces transformations, derived elsewhere, from the CS coupling 
scheme to the FS coupling scheme, or vice versa, related to the orbital symmetries [6]o 
and [42] o, appearing in the basis vectors (H). 

For the [6]o symmetry one obviously has : 

[6]o[33] F5 = [6]o[222] cs (Bl) 

either for IS=01 or 10. 

For the [42] o symmetry, sector IS=01, the Table |V| reproduces Table 7 of Ref. [^TJ with 



a phase change in columns 3 and 5, required by consistency with Ref. [26]. 
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In this Table, the column headings are 



ri 






= [42]o[321] cs 




= [42] [3111] os 




= [42] [222] cs 




= [42]o[21111]cs 



For the [42] q symmetry, sector IS=10, we reproduce in Table [V] the corresponding Table 
from Ref. |32|] by interchanging rows with columns and then reorder the rows. In this case, 
the notation is 



Vf 5 = 


[42]o[411] C5 


= 


[42] [33] C5 


^ s = 


[42] [2211]^ 


^ s = 


[42] [2211]^ 


$™ = 


[42] [l 6 ]cs 



The upper indices 1 and 2 take into account the two representations [2211] appearing in the 
inner product [222] x [33]. 



APPENDIX C: 

The calculation of the matrix elements of the Hamiltonian (|i~4"l) is based on the fractional 
parentage (cfp) technique described in Ref. fH| . For details, see also Ref. [jyj, chapter 10. 



In dealing with n particles the matrix elements of a symmetric two-body operator between 
totally (symmetric or) antisymmetric states ip n and ip' n reads 

77,(77, — 1 ) 

< V«l E >= 9 < V>n|K-l,n|< > (CI) 

i<j 1 

The matrix elements of V n -i tTl are calculated by expanding ip n and ip' n in terms of products 
of antisymmetric states of the first n-2 particles ip n -2 and of the last pair 02 
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^ = E P «WV- 2 («)0 2 (/?) (C2) 

a/3 

with a, (3 denoting the possible structures of and 02 and P a/ 3 the products of cfp 

coefficients in the orbital, spin-flavour and colour space states. In practical calculations, 



the colour space cfp coefficients are not required. The orbital cfp are taken from Ref. [|33 
Tables 1 and 2 by using the replacement rH 2 — ► s 4 p 2 and r 5 l —>■ s 5 p. The trivial ones are 
equal to one. The flavour-spin cfp for IS=01 are identical to the ^-matrices of Table 1 
Ref. |H| with [42]s[33]f in the column headings. For IS=10 they are the same as for IS=01 
but the column headings is [42] p [33] 5 instead of [42] 5 [33] p as above, and this is due to the 
commutativity of inner products of S n (see for example Ref. ||19||). The cfp used in the OC 



coupling are from Ref. @ Table 3, for [42] x [222] c -> [3111]oc and Table 5 of Ref. |3T 
for [42] x [222} c -> [222] oc and [42] G x [222] c -> [21111] 00 . 

In this way, after decoupling all degrees of freedom one can integrate out in the colour, 
spin and flavour space. The net outcome of this algebra is that any six-body matrix element 
becomes a linear combination of two-body orbital matrix elements, < V n >, < V v > and 

< Vtf >. The coefficients of < > are the same for IS=01 and 10, but the coefficients 
of < V v > are usually different. In both cases the coefficients of < V n < > are two times 
those of < V v >. We found that the two-body GBE matrix elements satisfy the relations 

< V n >~< V v > and < V v ' >— 2 < V n >. As an example, in Table [VI] we show the 
matrix elements obtained for IS=01. Except for < > and < s2s\V\pp > they are 
all negative, i.e. carry the sign of Eq. ([I]). 

In a harmonic oscillator basis the confinement potential matrix elements can be per- 
formed analytically. As an illustration, in Appendix D, we reproduce the results for all 
configurations required in these calculations. 

Finally, the kinetic energy matrix elements can be calculated as above, by writting the 
relative kinetic energy operator as a two-body operator 

^ = E||-^I>i) 2 = I>, (C3) 

with 
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T « = ikti + ®-L** (C4) 

Alternatively we can use an universal formula for the kinetic energy of harmonic oscillator 

states 

K.E. = ~[N + ~(n - l)]hcu + ~hcu (C5) 

where N is the number of quanta and n the number of particles. The last term is the kinetic 
energy of the center of mass. 



APPENDIX D: 



We work with the following single particle harmonic oscillator states : 

> =7r- 3 / 4 /r 3/2 exp(-r72/3 2 ) (Dl) 



|2s > = 2 1 / 2 3- 1 / 2 7T- 3 / 4 /T 3 / 2 (- - r —) exp (-r 2 /2/? 2 ) (D3) 



\p > rn = 8 1 / 2 3- 1 / 2 vr- 1 / 4 /?- 5/2 r exp (-r 2 /2/3 2 ) Y lm (D2) 

,3 - 2 
2 p 2 

In this basis the two-body matrix elements of the confining potential V c = Cr of Eq.flTo]) 

are 

< ss\V c \ss > = \ ~2CP (D4) 

V 7T 

< sp\V c \sp > = J- 1 — (D5) 

V 7T 3 

<sp\V c \ps> =-]fl^Y (D6) 

< ss\V c \(pp) L=0 > = -V3 < sp\V c \ps > (D7) 

< (pp)l=o\V c \(pp) l=0 > = ^ (D8) 

< s2s\V c \s2s > = \[^^- (D9) 

V 7T 12 

< ss\V c \s2s > = -4 hj-CP (D10) 



3tt 
1_C 
V^~2 



< s2s\V c \(pp) L=0 > = -±=~- (Dll) 
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TABLES 

TABLE I. Expectation values of the operators denned by Eqs. (|l|) and (|i~2] ) for all compatible 



symmetries [/]o[/]fs in the IS=(01) and (10) sector. < V x > is in units of C x and < V cm > in 
units of C cm - 







1=0 S=l 




1=1 S=0 


[fUfWs 


<Vx> 


< v cm > 


<Vx> 




[6] „ [33] fs 


-28/3 


8/3 


-8 


8 


[42] [3% s 


-28/3 


-26/9 


-8 


-4/3 


[42] [51] FS 


-100/3 


16/9 


-32 


16/9 


[42] [411] FS 


-28/3 


20/9 


-8 


44/9 


[42] [321] FS 


8/3 


-164/45 


4 


232/45 


[42] D [2211] FS 


68/3 


-62/15 


60 


42/5 
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TABLE II. Results of the diagonalization of the Hamiltonian (p3|) for IS=(01). Column 1 - 
the basis states, column 2 - diagonal matrix elements (GeV), column 3 - eigenvalues (GeV) for a 
5x5 matrix, column 4 - components of the lowest state. The results are for (3 = 0.437 fm. In 
columns 2 and 3, the quantity 2uin = 1.939 GeV is subtracted. 



state 


1 x 1 


5x5 


lowest state 
components 


|s 6 [6] [33] FS > 


2.346 


0.830 


-0.14031 


|sV[42] G [3% 5 > 


2.824 


1.323 


0.07747 


|sV [42] [51] FS > 


0.942 


2.693 


-0.96476 


|sV[42] [411] F5 > 


2.949 


3.049 


0.20063 


\(^/5/6s 5 2s - yi76 S y [6]o[3%s > 


3.011 


4.169 


0.05747 
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TABLE III. Same as Table 2 but for IS=(10). 



state 


1 x 1 


5x5 


lowest state 
components 


|s 6 f6l [33l F s > 


2.990 


1.356 


-0.12195 


\s 4 p 2 \42] \33] FS > 


3.326 


1.895 


0.08825 


|sV[42] [51] F5 > 


1.486 


3.178 


-0.96345 


|aV[42] [411] FS > 


3.543 


3.652 


-0.21644 


\(^5/6s 5 2s - v/IT^y [6] [3%s > 


3.513 


4.777 


0.04756 



TABLE IV. The unitary transformation between the CS and FS basis vectors of orbital sym- 
metry [42] o, isospin 1=0, and spin S=l. 



CS 



[42] [3% 5 
[42] [51] F5 
[42] [411] FS 
[42] [2211] FS 

[42] [321] F5 



9t/5 
36 



jh/E 
45 



9V10 
180 



11 
20 



18 
"45 



8VE 

36 



8^5 
45 



8V10 
180 



20 



29 
"45 



5\/2 
36 



5_\/2 
45 



170 
180 



20 



2V10 
45 



11 

36 



25 
"45 



25^ 
180 



5_\/5 
20 



lOy/5 
45 



20 

36 



25 
'45 



20^ 
180 



20 



8\/5 
45 
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TABLE V. Same as Table TV] but for S=0, 1=1. 



cs 



&8 



[42] [33] FS 
[42] [51] FS 
[42] [411] FS 
[42] [2211] FS 
[42] [321] F5 



'16 
15 
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TABLE VI. All one-meson exchange two-body matrix elements (GeV) for the sector IS=01 
evaluated at (3 = 0.437 fm. The remaining one is < sp\V F \ps >= — < ss\Vf\{pp)l=o > /V3. 



Two-body matrix 
elements 


F = TT 


F = i] 


F = rf 


< ss \Vf \ss > 


-0.108357 


-0.104520 


-0.189153 


< ss\V f \(pp)l=o > 


0.043762 


0.042597 


0.076173 


< sp \Vf \sp > 


-0.083091 


-0.079926 


-0.145175 


< (PP) L=0 1 V F | (jpp) L =0 > 


-0.081160 


-0.078594 


-0.142205 


< s2s\Vf\s2s > 


-0.069492 


-0.066963 


-0.121701 


< ss\Vf\s2s > 


-0.030945 


-0.030121 


-0.053863 


< s2s\V F \(pp) L =o > 


0.033309 


0.031753 


0.057499 



Figure Captions 
Fig. 1 

Projections of the lowest eigenvector in Table II onto AW and AA channels (in arbitrary 
units). 
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